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Abstract 

In this paper we study the effect of subtle changes in initial conditions on the evolu- 
tion of global quantities in two-dimensional Magnetohydrodynamic (MHD) turbulence. 
We find that a change in the initial phases of complex Fourier modes of the Elsasser 
variables, while keeping the initial values of total energy, cross helicity and Alfven ratio 
unchanged, has a significant effect on the evolution of cross helicity. On the contrary, 
the total energy and Alfven ratio are insensitive to the initial phases. Our simulations 
are based on direct numerical simulation using the pseudo-spectral method. 
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In fluid turbulence, the evolution of the velocity, u(x), at a given position or of a given 
Fourier component, u(k), is known to be sensitive to the details of the initial conditions, e.g. 
phases of u(k) (to be defined rigorously later). However, the evolution of global quantities 
like total energy are generally presumed to depend only on initial values of total energy 
and energy spectrum. The averaging over many modes appear to wash out the effects 
of the initial phases after a reasonably long time. The total energy in simulations with 
the same initial energy and spectrum, but with the modes chosen randomly, evolve along 
nearby trajectories; This is demonstrated in Fig. 1. It has been a common belief that in 
MHD turbulence also the evolution of global quantities, e.g. total energy and cross helicity, 
depends only on the initial values of global quantities and their spectra. The evolution 
of global quantities have earlier been studied by Ting et al. [|IJ, Matthaeus et al. 0, |3|], 
Biskamp and welter ||, Pouquet et al. |5[] in which they found dynamic alignment and 
various other phenomena. In this paper we show numerically that under certain conditions 
in two-dimensional (2D) magnetohydrodynamic (MHD) turbulence, the evolution of global 
quantities may not depend solely only on their initial values, but may depend significantly 
on more subtle features like the phases of complex Fourier modes (to be defined below) of 
the initially prescribed fields of the dynamical variables. In other words, a knowledge of the 
gross initial features as specified by the global quantities is not sufficient under all conditions 
to determine the evolution of global quantities. We only choose phases as a convenient way 
of demonstrating the inadequacy of specifying the initial global quantities and spectra alone 
for computing the evolution of certain global quantities. We find that the evolution of cross 
helicity shows a sensitivity to the initial phases in simulations with small values of initial 
cross helicity. 

The primary variables in MHD turbulence are the velocity field u and the magnetic field 
b. In our simulations we take mean magnetic field to be zero. We use Elsasser variables 



z ± = u ± b in our simulations. These variables denote fluctuations with 'positive' and 



'negative' velocity-magnetic field correlations. The relevant quadratic quantities are 
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the magnetic energy, 
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and the mean square magnetic vector potential 
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where ip is the magnetic vector potential, z + and z~ are complex conjugates of z + and z~ 
respectively. The total energy is E = (E + +E~)/2 and the cross helicity is H c = (E + —E~)/2. 
There are two well known dimensionless parameters: the normalised cross helicity a c = H c / E 
and the Alfven ratio ta = E u j The total energy, cross helicity, and the mean square vector 
potential are the three global inviscid invariants of the 2D MHD equation. 



We denote the complex Fourier modes z ± (k) by |z ± (k)| exp(i^), where 9^ are the phases 
of the modes. All the three global inviscid invariants E, H c , and A are independent of phases, 
while the Alfven ratio ta depends on the phase difference 0^ — 9^. Ting et al. [T|] found that 
the Alfven ratio affects the evolution of global quantities; it follows from their observations 
that the initial phase difference 9£ — 8^ would affect the global evolution. In this paper 
we demonstrate numerically that even keeping this initial phase difference fixed, change of 
absolute value of the initial phases 0£ affects the global evolution. 

In our simulations we investigate the effects of the initial phases on the subsequent total 
energy E, normalised cross helicity a c , and Alfven ratio r^. The temporal evolution of o c 
has been the subject of investigation in a number of earlier studies |], |], [3|, f|, 0. In several 
of these studies, a c has been observed to increase with time [0, 0, [|, |5J , a behaviour termed 
as dynamic alignment. However, Biskamp and Welter || observed in their simulations that 
the tendency towards dynamic alignment decreases with the increase in Reynolds number, 
and a c could even decrease at high enough Reynolds number J| . Ting et al. |lj too observed 
a few cases of decreasing a c for small values of initial a c and Ej A. In these earlier studies 
the effects of absolute phases have not been studied. 

We solve the 2-D incompressible MHD equations with hyperviscosity. The equations 
written in terms of the Elsasser variables, z + and z~ are 

^- =F (B . V) z ± + (z T .V) z ± = -Vp + i/±V 2 z ± + z/ T V 2 z T 

+ (u ± /k%) V 4 z± + (u T /k 2 eq ) V 4 zT (6) 

where u + and z/_ are related to the fluid viscosity (z/) and magnetic diffusivity (fi) by the 
relationship u ± = 1/2 [y ± fi). The last two terms in Eq. (6) include hyperviscosity v ± jk\ q 
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to damp out the energy at very high wave numbers. We choose v = \x = 5 x 10~ 4 for runs 
on a grid of size 512 x 512 and v = fi = 1CT 3 for runs on a grid of size 256 x 256. The 
hyperviscosity related parameter k eq is chosen to be 20 for runs on both the grids. The 
time step dt used for these runs is 10 -3 . The simulations are carried up to the final time 

t final = 50. 

We use the pseudo-spectral method [|, ^ § to solve the above equations in a periodic box 
of size 2tt x 2tt. In order to remove the aliasing errors arising in the pseudo-spectral method 
a square truncation is performed wherein all modes with \k x \ > N/3 or \k y \ > iV/3 are set 
equal to zero. The equations are time advanced using the second order Adam-Bashforth 
scheme for the convective terms and Crank-Nicholson for the viscous terms. In order to 
validate our code we used a simulation result of Pouquet et al. for comparison (see Fig. 
2). 

The simulations are performed for various initial sets of a c and ta values. The initial 
conditions are generated by first fixing ta, &c, E + , and E~ . The chosen value of ta determines 
the phase difference 9^ — 9^. The initial E and a c determine |z ± (k)|. Note that the absolute 
phase 9£ is still a free parameter. Only modes within the annular region 1/2 < |k| < 3/2 
are non-zero and each of the modes within this region receives equal energy (i.e., |z ± (k)| 2 = 
E /M, where M is the number of modes in the shell). The initial states are generated thus 
for only half the modes - the remaining half are conjugate to them. 

The phase sensitivity of the evolution of a c , E, and va are studied by comparison of pairs 
of simulations in which initial 9£ are different. We change the initial phases in two ways. In 
one case we change 9^ uniformly for all the modes by an amount A, while in the other case 
the phases are changed by using a different random seed in the random number generator. 
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The initial global quantities E, H c , ta, and their spectra remain unchanged under these 
phase changes. The evolution of <j c for a variety of initial a c , ta and A values are shown 
in Table 1 (pairs of simulations are shown together; for example, mhdl differs from mhdl* 
only in that its initial z ± (k) fields have to be shifted from the latter's by A = 0.4). 

The N = 512 runs with tfi na i = 50 are very time intensive. Hence only the small initial 
a c runs, which we found to be sensitive to the phases, were carried out for N=512. A large 
number of runs were performed on N=64 to explore a wider range of initial conditions. All 
these results showed behaviour consistent with the results discussed below which are based 
on the high resolution runs N=256 and 512. 

In our simulations, the small a c runs showed the most significant dependence on initial 
phase 9£. For a c = 0.1 and ta = 1.5, we choose A to be 0.0 and 0.4 in mhdl and mhdl* 
respectively. It is seen in Fig. 3 that phase shifting has a marked effect on the evolution of 
a c . For A = 0.0 (mhdl) a c increases from its initial value of 0.1 to its final value of 0.129, 
whereas for A = 0.4 (mhdl*), a c decreases to a final value of 0.06. The total energy (Fig. 
4) and the Alfven ratio (Fig. 5) do not appear to be affected much by the phase shift. We 
also compare two simulations (mhdl and mhdl** in Table 1.) in which the initial phases are 
generated using different random seeds. For these cases also the effect on the evolution of 
a c (Fig. 3) is significant, but the corresponding effects on the evolution of the total energy 
(Fig. 4) and ta (Fig. 5) are not noticeable. We also studied the effects of initial phases for 
the same initial a c , but with a large initial 7-^(5.0). The runs mhd2 and mhd2* show the 
effects of changing A, while mhd2 and mhd2** show the effects of different random number 
generator seeds. The results obtained for this case are similar to the run for initial condition 
with ta = 15. In Fig. 6 it is seen that for initial value of A = 0.0 (mhd2), a c increases 
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and for initial A = 0.3, a c decreases. The effect of changing A on the total energy (Fig. 7) 
and ta (Fig. 8) is seen to be small. Similar results are obtained if we change the seed of the 
random number generator (compare mhd2 and mhd2**). Hence, a c (Fig. 6) is sensitive to 
the change in the initial phases whereas the total energy (Fig. 7) and ta (Fig. 8) are not 
sensitive. Earlier, Ting et al. [0] had observed a decrease in a c for small initial a c . 

We also perform runs at higher initial values of a c (mhd3 and mhd3* in Table 1.). The 
effect of phase shifting for initial a c = 0.5 and ta = 5.0 is shown in Figures 9,10,11. It is seen 
in Fig. 9 that the changes in evolution caused by phase shifting are relatively smaller for 
high initial a c values as compared to small initial a c discussed above. From Fig. 10 it can be 
seen that the effect of A on total energy is also small and ta (Fig. 11) remains insensitive to 
the change in A. We have performed more runs than have been shown here and in all cases 
cr c was seen to increase for high a c values. This increase in a c is consistent with simulations 
performed earlier for high a c values ]T|, |3|, f|, [5| . 

From the numerical results presented here we conclude that phases of the initial modes 
play an important role in the evolution of a c , atleast for cases with small initial a c values. 
For higher values of <r c , phases do not appear to affect the evolution of a c by any significant 
amount. In all the runs the total energy and ta were seen not to have any significant 
dependence on the phases. 

The origin of the phase effects discussed here is not clear at this moment. We need to 
examine the evolution more carefully before reaching any definite conclusion. These studies 
could find applications in understanding the solar wind observations in which a c has been 
observed to decrease ||, [TO], [H], 

It has been demonstrated in the paper that the evolution of normalised cross helicity 

7 



is significantly affected by subtle features of the initial condition especially at low initial 
cross helicities. This observation will require us to be more circumspect in drawing conclu- 
sions based on arbitrary initial conditions and to exercise more care in choosing the initial 
conditions in MHD turbulence. 
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FIGURE CAPTIONS 



Figure 1. Energy evolution for fluids. The two cases shown here differ only in the phases of 
the initial conditions. The modes in the runs were |u(k)| exp^ 6>k+A ) with A = 0.0 and 0.4. 

Figure 2. Time evolution of kinetic (solid line) and magnetic (dashed line) energies for the 
O-T vortex using a set of parameters from Pouquet et al. (Pouquet et al's results are shown 
by diamonds and crosses). Here v = 2.5 x 10 -3 and hyperviscosity is zero. 

Figure 3. Normalised cross helicity (a c ) evolution for initial a c = 0.1, ta = 1.5. The curves 
shown correspond to mhdl, mhdl*, and mhdl** in Table 1. 

Figure 4. Evolution of total energy (E) for same initial conditions as in Fig. 3. See Table 1. 
for description of mhdl, mhdl*, and mhdl**. 

Figure 5. Evolution of Alfven ratio (r A ) for same initial conditions as in Fig. 3. Look at 
Table 1. for description of mhdl, mhdl*, and mhdl**. 

Figure 6. Normalised cross helicity (a c ) evolution for initial a c = 0.1, ta = 5.0. The curves 
shown correspond to mhd2, mhd2*, and mhd2** in Table 1. 
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Figure 7. Evolution of total energy (E) for same initial conditions as in Fig. 6. Look at 
Table 1. for description of mhd2, mhd2*, and mhd2**. 

Figure 8. Evolution of Alfven ratio (r A ) for same initial conditions as in Fig. 6. Look at 
Table 1. for description of mhd2, mhd2*, and mhd2**. 

Figure 9. Normalised cross helicity (<r c ) evolution for initial a c = 0.5 and = 5.0. The 
curves shown correspond to mhd3, mhd3* in Table 1. 

Figure 10. Evolution of total energy (E) for same initial conditions as in Fig. 9. Look at 
Table I. for description of mhd3, mhd3*. 

Figure 11. Evolution of Alfven ratio for same initial conditions as in Fig. 9. Look at Table 
1. for description of mhd3, mhd3*. 
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Table 1: Initial values of the random number generator seed A, a c and for runs performed 
on grid of size N x N. The initial and the final values (at tfi na i = 50) of a c are also shown. 
Run N Seed ta A a c (t = 0) a c (t = 50) a c increases/decreases 



mhdl 512 50 1.5 0. 0.1 



mhdl* 512 50 1.5 0.4 0.1 



mhdl** 512 575 1.5 0. 0.1 



mhd2 512 50 5.0 0. 0.1 



mhd2* 512 50 5.0 0.3 0.1 



mhd2** 512 575 5.0 0. 0.1 



mhd3 256 50 5.0 0. 0.5 



mhd3* 256 50 5.0 0.4 0.5 



mhd4 64 50 1.5 0. 0.1 



mhd4* 64 50 1.5 0.6 0.1 



mhd4** 64 50 1.5 1.0 0.1 



mhd5 64 50 2.0 0. 0.1 



mhd5* 64 50 2.0 0.2 0.1 



mhd5** 64 50 2.0 0.4 0.1 



mhd6 64 50 1.0 0. 0.5 



mhd6* 64 50 1.0 0.2 0.5 



mhd6** 64 50 1.0 0.4 0.5 



0.06 
0.13 
0.20 
0.22 
0.05 
0.13 
0.87 
0.88 
0.02 
0.34 
0.17 

- 0.02 

- 0.01 
0.22 
0.79 
0.74 
0.72 



decreases 
increases 
increases 
increases 
decreases 
increases 
increases 
increases 
decreases 
increases 
increases 
decreases 
decreases 
increases 
increases 
increases 
increases 
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Fig. 3 



15 




Fig. 4 
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Fig. 9 
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Fig. 10 
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